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We investigate the relationship between the blood flow velocities (BFV) in the middle cerebral 
arteries and beat-to-beat blood pressure (BP) recorded from a finger in healthy and post-stroke 
subjects during the quasi-steady state after perturbation for four different physiologic conditions: 
supine rest, head-up tilt, hyperventilation and CO2 rebreathing in upright position. To evaluate 
whether instantaneous BP changes in the steady state are coupled with instantaneous changes in the 
BFV, we compare dynamical patterns in the instantaneous phases of these signals, obtained from 
the Hilbert transform, as a function of time. We find that in post-stroke subjects the instantaneous 
phase increments of BP and BFV exhibit well pronounced patterns that remain stable in time for all 
four physiologic conditions, while in healthy subjects these patterns are different, less pronounced 
and more variable. We propose a new approach based on the cross-correlation of the instantaneous 
phase increments to quantify the coupling between BP and BFV signals. We find that the maximum 
correlation strength is different for the two groups and for the different conditions. For healthy 
subjects the amplitude of the cross-correlation between the instantaneous phase increments of BP 
and BFV is small and attenuates within 3-5 heartbeats. In contrast, for post-stroke subjects, this 
amplitude is significantly larger and cross-correlations persist up to 20 heartbeats. Further, we 
show that the instantaneous phase increments of BP and BFV are cross-correlated even within a 
single heartbeat cycle. We compare the results of our approach with three complementary methods: 
direct BP-BFV cross-correlation, transfer function analysis and phase synchronization analysis. Our 
findings provide new insight into the mechanism of cerebral vascular control in healthy subjects, 
suggesting that this control mechanism may involve rapid adjustments (within a heartbeat) of the 
cerebral vessels, so that BFV remains steady in response to changes in peripheral BP. 

PACS numbers: 05.40.-a,87.19.Hh,87.10.+e 



I. INTRODUCTION 



Cerebral autoregulation (CA) is the ability of cerebral 
blood vessels to maintain steady cerebral perfusion in re- 
sponse to fluctuations of systemic blood pressure (BP), 
postural changes or metabolic demands. This regula- 
tory mechanism is known to operate over a range of 
blood pressure values (e.g. 80 - 150 mm Hg) and 
on time scales of a few heartbeats 0- The long-term 
CA compensates for chronic BP elevations and metabolic 
demands Ischemic stroke is associated with an im- 
pairment of autoregulation |^ |^ , that may permanently 
affect cerebrovascular reactivity to chemical and blood 
pressure stimuli HQ. With impaired CA, the cerebral 
blood flow follows BP fluctuations, posing a risk of in- 
sufficient blood supply to the brain during transient re- 
ductions in peripheral BP. Therefore, evaluating the ef- 
fectiveness of CA is of great interest, given the clinical 
implications. 



*VN and PChI contributed equally to this work. VN designed 
the clinical study, provided data and guidance of the physiological 
aspects. PChI proposed the instantaneous phase cross-correlation 
method, and guided the computational aspects. 



Traditional experimental methods evaluating the 
mechanism of CA require time-consuming invasive proce- 
dures and are focused on long term BP and blood 
flow velocities (BFV) characteristics such as the mean, 
lacking descriptors of temporal BP-BFV relationship. To 
address this problem, an alternative "dynamic" approach 
has been proposed .10] to quantify CA using transcra- 
nial Doppler ultrasound during the transient responses 
in cerebral BFV to the rapid BP changes induced ex- 
perimentally by thigh cuff inflation, Valsalva maneuver, 
tilt-up, or a change in posture [sLlllf. The autoregulation 
indices derived from this approach may be more sensitive 
to indicators of hypoperfusion after stroke ^3 • 

The analytic methods evaluating the dynamics of cere- 
bral autoregulation are currently based on mathemati- 
cal modeling and Fourier transform analysis [l^ . The 
Fourier transform based transfer function method has 
been widely used 0. This method estimates the rela- 
tive cross-spectrum between BP and BFV signals in the 
frequency domain. Dynamic indices of autoregulation, 
based on Fourier transform methods, presume (i) signal 
stationarity (i.e., the mean and standard deviation of the 
signal are stable and remain invariant under a time shift) 
and (ii) a linear BP-BFV relationship. However, physio- 
logic signals are often nonstationary reflecting transient 
responses to physiologic stimuli Till. The effect of this 
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nonstationarity on the results obtained from the trans- 
fer function analysis has not been carefully assessed in 
previous studies. 

Here we investigate the dynamics of the BP-BFV rela- 
tionship when the system reaches the quasi-steady state 
after an initial perturbation. While studies traditionally 
have focused on the response in BFV to transient changes 
in BP 3], we hypothesize that spontaneous physiologic 
fluctuations during the quasi-steady state, which is char- 
acterized by absence of physiologic stimuli or constant 
level of stimulation, may also contain important infor- 
mation about the CA mechanism. 

Our focus on fluctuations in the BP and BFV is mo- 
tivated by previous work which has demonstrated that 
physiologic fluctuations contain important information 
about the underlying mechanisms of physiologic control. 
Robust temporal organization was reported for the fluc- 
tuations characterizing cardiac dynamics (inter-beat in- 
tervals) p, El El El El mjMJlpLresraratory dy- 
namics (inter-breath intervals) 123. |24L 25L l2a. |27[ . loco- 
motion (gait, fore- arm motion) |28l. l2£ , ISOl l3lL 132. l33l|. 
and brain dynamics 0, 1^1 Hal ■ Moreover, it has been 
demonstrated that physiologic fluctuations carry infor- 
mation reflecting the coupling between different physio- 
logic systems, e.g., correlations in the heartbeat change 
with physical activity [sl HI] , with wake and sleep [33 , 
during different sleep sta ges |40l liH li^ . li^ . and even 
different circadian phases |4J|. BP and BFV signals are 
impacted by the heartbeat, and thus one can expect that 
fluctuations in BP and BFV may reflect modulation in 
the underlying mechanisms of control. Previous stud- 
ies have focused on the transitional changes in BP and 
BFV in response to abrupt perturbation of the physio- 
logic state, e.g., rapid switch from supine to tilt. In con- 
trast, we focus on the dynamical characteristics of BP 
and BFV signals after the initial perturbation, when the 
system has reached a quasi-steady state during which 
there is no change in physiologic stimuli. Further, we 
hypothesize that certain dynamical characteristics of the 
fluctuations in BP and BFV at the steady state, and how 
these characteristics change for different physiologic con- 
ditions, may reflect aspects of the underlying mechanism 
of CA. For example, under normal CA the fluctuations 
of the BFV in healthy subjects at the steady state may 
relate to high frequency adjustments (even within a sin- 
gle heartbeat) of the diameter of the cerebral blood ves- 
sels, while loss of CA after stroke may lead to impaired 
vascular dilation or contraction associated with reduced 
fluctuations in BFV. 

To test this hypothesis, we measure BP and BFV sig- 
nals from healthy and post-stroke subjects during four 
physiologic conditions: supine, tilt, hyperventilation, and 
CO2 rebreathing in upright position. We apply several 
complementary methods to quantify the dynamical BP- 
BFV relationship in these quasi-steady conditions: trans- 
fer function analysis, cross-correlation and phase syn- 
chronization analyses, and we compare these methods 
with a new approach of cross-correlation between the in- 



stantaneous phase increments of BP and BFV signals. 
Interactions between peripheral circulation (beat-to-beat 
BP) and cerebral vasoregulation [BFV in the middle cere- 
bral artery (MCA)] can be modeled as dynamic synchro- 
nization of two coupled nonlinear systems. Specifically, 
We hypothesize that the CA mechanism may also involve 
adjustments in the cerebral vascular tone to spontaneous 
changes in BP that may be present within a single heart- 
beat even when the system is in the steady state and 
there are no significant changes in the mean blood pres- 
sure. 



The synchronization phenomenon was first observed 
by Hugenii for two coupled pendulum clocks [45l |. and 
since then it has been found in many physical and bio- 
logical systems where two or more coupled subsystems 
interact jH EHil El HI HI HI . Alternatively, the 
synchronization may also be triggered by the influence 
of external noisy or regular fields [sl |3- I^i recent 
years, the concept of synchronization has been widely 
used to study the coupling of oscillating systems, leading 
to the discovery of phase synchronization in non-identical 
coupled systems in which the instantaneous phases are 
synchronized, while their instantaneous amplitudes re- 
main uncorrelated [SI IH Ell. Such phase syn- 
chronization has been empirically discovered in a range 
of physical and physiological systems [sl |63| . Specifi- 
cally, studies have found coupling between the cardiac 
rhythm and other systems: phase synchronization was 
observed between the heartbeat and respiration during 
normal conditions [sij and during respiratory sinus ar- 
rhythmia '62', 'gs'I , change of cardiorespiratory phase syn- 
chronization with new-borns' age .64] and with the level 
of anaesthesia [b^. Phase analysis methods have been 
used to probe spatial synchronization of oscillations in 
blood distributio n sy stems , between cortical centers 
during migraine |67| . as well as between certain brain 
areas and muscle activity of the limbs . 



In this study we evaluate the time-domain characteris- 
tics of both the amplitudes and the instantaneous phases 
of the BP and BFV which can be considered as two in- 
teracting subsystems within the CA mechanism. To de- 
termine the characteristics of the coupling between BP 
and BFV in healthy subjects and how they change with 
stroke, we analyze the cross-correlation between the in- 
stantaneous phase increments of these two signals. We 
find that this cross-correlation is much stronger for post- 
stroke subjects, indicating increased synchronization be- 
tween BP and BFV, which suggests impaired mechanism 
of the CA. We compare the results of the instantaneous 
phase increment cross-correlation analysis with those ob- 
tained from several complementary methods including 
the transfer function, cross-correlation and phase syn- 
chronization analyses. 
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Variable Dcmogiapliic characteristics 



Control Stroke 
Mon/Womon 4/7 7/6 

Age(mean±SD) 48.2 ± 8.7 52.8 ± 7.1 

Race W/AA 10/1 12/1 

Stroke side — 5/8 

(Right/Left) 

Group Supine Tilt Hyperventilation CO2 rcbreathing Statistics 

(mcaniSD) Control Stroke Control Stroke Control Stroke Control Stroke (P values) 

BP (mm Hg) 96.4 ± 20.9 101.0 ± 20.6 93.1 ± 20.1 105.9 ± 22.5 94.2 ± 19.7 105.4 ± 22.0 97.6 ± 21.7 108.8 ± 22.7 0.66* 0.0005 + 
BFV-MCAR/ 66.0 ± 18.7 55.2 ± 18.0 57.3 ±18.1 49.7 ± 18.5 40.3 ± 15.7 39.8 ± 14.9 54.7 ± 21.9 56.3 ± 21.5 < 0.0001* 0.77* 
Normal side(cm/s) 

BFV-MCAL/ 63.5 ±19.6 51.1 ± 19.0 54.8 ±17.1 51.5 ± 19.3 40.8 ± 14.9 41.0 ± 17.9 54.6 ± 21.5 57.4 ± 23.3 < 0.0001* 0.008+ 
Stroke side(cm/s) 

C02(mm Hg) 33.5 ± 6.0 37.7 ± 4.9 32.0 ± 3.6 32.5 ± 2.5 21.0 ± 4.5 23.5 ± 8.0 34.6 ± 7.2 33.2 ± 6.5 < 0.0001* 0.28 + 
CVR**-MCAR 1.54 ± 0.45 1.96 ± 0.54 1.75 ± 0.53 2.28 ± 0.93 2.68 ± 1.10 3.00 ± 1.28 1.97 ± 0.61 2.17 ± 0.83 0.0001* 0.007+ 
/Normal side 

CVR-MCAL 1.59 ± 0.56 1.97 ± 0.65 1.79 ± 0.61 2.29 ± 0.83 2.57 ± 0.90 3.04 ± 1.60 1.96 ± 0.56 2.19 ± 0.94 0.0001* 0.01 + 
/Stroke side 

* P value between physiologic conditions comparisons 
f P value between groups comparisons 
** CVR (cerebral vascular resistance) is defined as mean BP/BFV 



TABLE I: Demographic characteristics. 



II. METHODS 



A. Study Groups 



We obtain data from the Autonomic Nervous System 
Laboratory at the Department of Neurology at The Ohio 
State University and from the SAFE (Syncope and Falls 
in the Elderly) Laboratory at the Beth Israel Deaconess 
Medical Center at Harvard Medical School. All subjects 
have signed informed consent, approved by the Institu- 
tional Review Boards. Demographic characteristics are 
summarized in Table Control group: 11 healthy sub- 
jects (age 48.2 ± 8.7 years). Stroke group: 13 subjects 
with a first minor ischemic stroke ( > 2 months after 
acute onset) (age 52.8 ± 7.1 years). Post-stroke subjects 
have a documented infarct affecting < 1/3 of the vascu- 
lar territory as determined by MRI or CT with a minor 
neurological deficit (Modified Rankin Score scale < 3). 
The side of the lesion is determined by neurological eval- 
uation and confirmed with MRI and CT. The lesion is in 
the right hemisphere in 5 of the subjects and in the left 
hemisphere in 8 of the subjects. Normal carotid Doppler 
ultrasound study is required for participation. Patients 
with hemorrhagic strokes, clinically important cardiac 
disease including major arrhythmias, diabetes and any 
other systemic illness are excluded. All subjects are care- 
fully screened with a medical history, physical and labo- 
ratory examination. 



B. Experimental Protocol 

All subjects have participated in the following experi- 
mental protocol: 

• Baseline supine rest — normal breathing (normocap- 
nia): subject rests in supine position for 5 minutes on a 
tilt table; 

• Head-up tilt — upright normocapnia: The tilt table is 
moved upright to an 80° angle. The subject remains in 
upright position for 5 minutes and is breathing sponta- 
neously; 

• Hyperventilation — upright hypocapnia: the subject is 
asked to breathe rapidly at «1 Hz frequency for 3 min- 
utes in an upright position. Hyperventilation induces 
hypocapnia (reduced carbon dioxide) , which is associated 
with cerebral vasoconstriction; 

• CO2 rebreathing — upright hypercapnia: The subject 
is asked to breath a mixture of air and 5% CO2 from 
rebreathing circuit at a comfortable frequency for 3 min- 
utes in an upright position. CO2 rebreathing increases 
carbon dioxide above normal levels and induces hyper- 
capnia, which is associated with vasodilatation. 

The mechanism of CA is at least partially related to 
the coupling between metabolic demands and oxygen 
supply to the brain Q]. Carbon dioxide (CO2) is one 
of the most potent chemical regulators of cerebral va- 
soreactivity. Head-up tilt provides both pressure and 
chemical stimulus — BFV and CO2 decline in upright 
position, reflecting the change in intracranial pressure 
and shifting autoregulatory curve towards lower BP val- 
ues. There is a linear relationship between CO2 values 
and cerebral blood flow: hypocapnia (through hyper- 
ventilation) causes vasoconstriction and thus decreases 
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FIG. 1: BP and BFV signals during CO2 re-breathing after a band-pass Fourier filter in the range [0.05Hz, fOHz] and normal- 
ization to unit standard deviation; (a-b) for a healthy subject, (c-d) for a post-stroke subject. 



the blood flow, hypercapnia (through CO2 rebreathing) 
causes vasodilatation and increases the blood flow in the 
brain 



C. Data Acquisition 

We perform experiments in the morning or more than 
2 hours after the last meal. We measure the electrocar- 
diogram from a modified standard lead II or III using 
a SpaceLab Monitor (SpaceLab Medical Inc., Issaquah, 
WA) . We record beat-to-beat BP from a finger with a Fi- 
napres device (Ohmeda Monitoring Systems, Englewood 
CO), which is based on a photoplethysmographic vol- 
ume clamp method. During the study protocol, we ver- 
ify BP by arterial tonometry. With finger position at 
the heart level and temperature kept constant, the Fi- 
napres device can reliably track intraarterial BP changes 
over prolonged periods of time. We measure the res- 
piratory waveforms with a nasal thermistor. We mea- 
sure CO2 from a mask using an infrared end tidal vol- 
ume CO2 monitor (Datex Ohmeda, Madison WI). We 
insonate the right and left MCAs from the temporal win- 
dows, by placing the 2-MHz probe in the temporal area 
above the zygomatic arch using a transcranial Doppler 
ultrasonography system (MultiDop X4, DWL Neuroscan 
Inc, Sterling, VA). Each probe is positioned to record 
the maximal BFV and is fixed at a desired angle us- 
ing a three-dimensional positioning system attached to 
the light- metal probe holder. Special attention is given 
to stabilizing the probes, since their steady position is 
crucial for reliable, continuous BFV recordings. BFV 
and all cardiovascular analog signals are continuously ac- 
quired on a beat-to-beat basis and stored for off-line post- 
processing. We visually inspect the data and remove oc- 
casional extrasystoles and outlier data points using lin- 



ear interpolation. We use the Fourier transform of the 
Doppler shift (the difference between the frequency of the 
emitted signal and the echo frequency of the reflected sig- 
nal) to calculate BFV. BFVs in the MCA correlate with 
invasive measurements of blood flow with xenon clear- 
ance J^ser Dopp ler flux and positron emission tomogra- 
phy lii, [73, EJ • Since the MCA diameter is relatively 
constant under physiological conditions jT^j , BFV can be 
used for blood flow estimates. 



D. Statistical Analyses 

We use the multiple analysis of variance (MANOVA) 
with 2x4 design for the two groups (control and stroke) 
and for four physiologic conditions (supine, tilt, hy- 
perventilation, and CO2 rebreathing) with subjects as 
nested random effects (JMP version 5 Software Analysis 
Package, SAS Institute, Gary, NC). For each group and 
condition, we calculate: (i) the mean BP, BFV, cerebral 
vascular resistance (CVR, calculated by mean BP/BFV) 
from the right and the left MCAs, and CO2 (see Ta- 
ble QJ; (ii) gain, phase and coherence from transfer func- 
tion analysis (see Table Hi)) : and (iii) the direct cross- 
correlation, the phase synchronization, and the cross- 
correlation of instantaneous phase increments of BP and 
BFV (see Table IHT)! . We analyze BFV on the stroke side 
in the right MCA in 5 patients and in the left MCA in 
8 patients. In our comparative tests, we consider the 
side opposite to the stroke side as the "normal" side, as 
we did not a priori know whether the side opposite to 
the stroke side would exhibit normal or perturbed be- 
havior. In the group comparison, we compare the stroke 
side BFV for the stroke group to the BFV in the left side 
MCA for the control group because the majority of post- 
stroke subjects had stroke on their left side. We note 
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FIG. 2: Cross-correlation function C(r) for the BP and BFV signals during four physiologic conditions: (a) for the same 
healthy subject and (b) for the same post-stroke subject as shown in Fig. BP and BFV signals are pre-processed using a 
band-pass Fourier filter in the range [0.05Hz, lOHz] and are normalized to unit standard deviation before the analysis. Since 
BFV precedes BP, the maximum value Cmax in the cross-correlation function C(t) is located not at zero lag but at r ~ —0.1 
sec. 



that our comparative tests between the stroke side BFV 
for post-stroke subjects and the right side BFV for the 
control group show similar results. We compare the nor- 
mal side (opposite to the stroke side) BFV for the stroke 
group to the BFV in the right side MCA for the control 
group. 

Method 1: Transfer Function Analysis 

We first normalize the BP and BFV signals to unit 
standard deviation to obtain the respective signals P{t) 
and V{t). We then calculate their respective Fourier 
transforms V{f) and P{f)- In the frequency domain, 
the coherence function 7^(/) is defined as 



\Spvif)? 



Spp{f)Svv{f)' 



(1) 



where Svvif) = \V{f)\\ Spp{f) = \P{f)\^ and 
Spvif) = P*{f)y{f) are the power spectra of V{t), 
P{t) and the cross-spectrum of V{t) and P{t), respec- 
tively. The value of the coherence function 7^(/) varies 
between and 1. The transfer function H{f) is defined 
as 



Spy if) 
Sppif)' 



(2) 



From the real part Hji{f) and imaginary part Hj{f) of 
the transfer function, we can obtain its amplitude (also 
called gain) 



\H{f)\^[HUf)+Hj{f)] 



1/2 



(3) 



We note that the phase <&(/) is a frequency domain char- 
acteristic of the cross-spectrum between two signals, and 
is different from the instantaneous phase in the time do- 
main we discuss in Method 3 and Method 4- 

Method 8: Cross-correlation Analysis 

To test dynamical aspects of the mechanism of CA, 
we investigate the cross-correlation between BP and 
BFV signals and how this cross-correlation changes un- 
der stroke. Although we consider segments of the BP 
and BFV recordings during the quasi-steady state, where 
we have constant level (or absence) of physiologic stim- 
uli, these signals may still exhibit certain trends in the 
mean value. So, to eliminate these trends, we first pre- 
process the BP and BFV signals applying a band-pass 
filter (lOHz > / > 0.05Hz) in the frequency domain. To 
be able to compare signals with different amplitude, we 
next normalize the band-passed signals to unit standard 
deviation (Fig.QJ. Finally, we perform a cross-correlation 
analysis estimating the cross-correlation function C(t) 
for a broad range of values for the time lag r, where 
C(t) is defined as: 



C(r) 



mt)^{P)){V{t + T)-{V))) 



(5) 



and its phase 



and (...) and cr denote the mean and standard deviation 
of a signal. 

Results of our cross-correlation analysis for one healthy 
subject and one post-stroke subject during all four phys- 
iologic conditions are shown in Fig. |21 and are discussed 

in sec. rrrrci 



$(/) EEarctan[i/,(/)/ffH(/)]. 



(4) Method 3: Phase Synchronization Analysis 
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FIG. 3: (a) The relative phases ipBFv{t) — tpBp{t) for heahhy subjects (a) and post-stroke subjects (b), where ip{t) is the 
instantaneous phase of a signal. The phase difference of BFV and BP for both groups may fluctuate around a constant value 
or jump between different constant values. The distributions of V&i,! = (27r)~^[</3sFv(t) — V'sp(t)] mod 1 for healthy subjects 
and for post-stroke subjects are shown in (c) and (d), respectively. The number of bins in the histogram is A'' = 40. 



We study beat-to-beat BP-BFV interaction during 
quasi-steady state conditions (supine rest, upright tilt, 
upright hyperventilation and CO2 rebreathing employing 
a phase synchronization method. We first apply a high- 
pass (/ > 0.05Hz) and a low-pass (/ < lOHz) Fourier 
filter to the BP and BFV signals. The high-pass filter is 
used to reduce nonstationarity related to slow trends in 
the mean of the signals. The low-pass filter is used to re- 
move high frequency random fluctuations in the signals. 
Next we perform the Hilbert transform which for a time 
series s (t) is defined as jH HI [zl IH [zl IH 

m ^ Ip r ^dr, (6) 

where P denotes Cauchy principal value. s{t) has an 
apparent physical meaning in Fourier space: for any pos- 
itive (negative) frequency /, the Fourier component of 



the Hilbert transform s{t) at this frequency / can be ob- 
tained from the Fourier component of the original signal 
s{t) at the same frequency / after a 90° clockwise (anti- 
clockwise) rotation in the complex plane. For example, 
if the original signal is sn\(ujt)^ its Hilbert transform will 
become cos{ujt). For any signal s(t) one can always con- 
struct its "analytic signal" S [sl IslFzl Fzl 1?^ . which 
is defined as 

S = s{t) + i~s{t) ^ A{t)e''^^'\ (7) 

where A[t) and Lp{t) are the instantaneous amplitude and 
instantaneous phase of s{t), respectively. Application of 
the analytic signal approach to heartbeat dynamics has 
been shown in |l5l l75l |. The instantaneous amplitude 
A{t) and the instantaneous phase Lp(t) are instantaneous 
characteristics of a time series s{t), and present different 
aspects of the signal. For a pure sinusoid, A{t) is a con- 
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FIG. 4: (a) The relative phases ipBFv{ti) — ipBp{t2) for a healthy subject and a post-stroke subject. We find that the relative 
phases between BFV and BP depend on the time difference t2 — ti between these two signals. Accordingly, the distributions of 
^1,1 = {2'K)~^[ipBFv{ti) — 'fiBp{t2)] mod 1, shown in (b) and (c) for a healthy subject and a post-stroke subject respectively, 
also depend on the time difference between BFV and BP. We note that, our choice of time difference t2 — t\ — 0.4 sec is 
arbitrary and does not carry any specific physiologic meaning. The number of bins in the histogram is TV = 40. 



stant and ip[t) increases linearly in time: the amplitude 
quantifies the strength of the oscillation and the slope of 
the straight line formed by the increasing phase quanti- 
fies how fast is the oscillation. For more complex signals, 
both A(t) and (fit) may display complicated forms. Fur- 
ther, we note that the instantaneous phase ip{t) is dif- 
ferent from the transfer function phase $(/): <p{t) is a 
time domain characteristic of a single signal, while $(/) 
is a cross-spectrum characteristic of two signals in the 
frequency domain. 

From the definition of the Fourier transform one can 
find that the mean of a given signal sit) is proportional 
to the Fourier component of s{t) at the frequency / = 0. 
After applying the high-pass filter (/ > 0.05 Hz), the 
Fourier component of s(i) at / = is filtered out, and 
correspondingly, the mean of the filtered signal becomes 
zero. Furthermore, from the definition of the Hilbert 
transform, one can find that the Hilbert transforms of 
two signals s{t) and s{t) + c (where c 7^ is a constant) 



are identical, although the two original signals have ob- 
viously different means due to the constant c. Since both 
the instantaneous phase ip{t) and amplitude A{t) depend 
on the original signal as well as on the Hilbert transform 
of the original signal (see Eq.[7I), the instantaneous phase 
and amplitude for s{t) and s(i) -I- c will be also different. 
However, when we let both signals s(t) and s{t) -\- c pass 
a high-pass filter (/ > 0.05 Hz), the mean of both signals 
becomes zero, and both signals will have identical phases 
and amplitudes. Thus, the phase and amplitude of a sig- 
nal does not depend on its mean, and is uniquely defined 
after the signal is processed with a high-pass filter. 

Following [53 we estimate the difference between 
the instantaneous phases of the BFV and BP signals 
ipBFv(t) — ipBp(t) at the same time t. We then obtain 
the probability density p for the quantity 

^-14 = {2TT)-'[ipBFv{t) - ifiBpit)] mod 1, (8) 

following the approach presented in |68j| . To quantify the 
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FIG. 5: Presentation of the BP and BFV signals vs. their Hilbert transforms (a-b) and their corresponding instantaneous phase 
increment A(/? during the CO2 re-breathing condition (c-d) for the same data from a healthy subject as shown in Fig. ^-b. BP 
and BFV signals vs. their Hilbert transforms (e-f) and their corresponding instantaneous phase increment A<p during the CO2 
re-breathing condition (g-h) for the same data from a post-stroke subject as shown in Fig.^-d. Repetitive temporal patterns 
associated with each heartbeat in Aif for the BP signal from a healthy subject (c) are not matched by corresponding patterns 
in the BFV signal (d), reflecting active cerebral vascular regulation. In contrast, periodic patterns in Aip of the BP signal 
from a post-stroke subject (g) are matched by practically identical patterns in A<p of the BFV signal (h), indicating dramatic 
impairment of cerebral vascular tone with higher vascular resistance after minor ischemic stroke. 



shape of the probability density p, which can be very 
broad (uniform) for non-synchronized signals and very 
narrow (peaked) for strongly synchronized signals, we 
utilize the index p j68|. The index p is based on the 
Shannon entropy and is defined as 

P = (-Smax - '5')/5'max, (9) 

where S = — X^fcLi Pfclnpfc is the Shannon entropy of the 
distribution of and 5'max = InA'^ corresponds to the 
uniform distribution, where N is the number of bins. In 
our calculations, we choose N = 40. For the normalized 
index p, we have < p < 1, where p = corresponds 
to a uniform distribution (no synchronization) and p = 1 



corresponds to a Dirac-like distribution (perfect synchro- 
nization). 

Results of our synchronization analysis for several 
healthy and post-stroke subjects are shown in Figs. |3| 
and^ and are discussed in Sec. IIII Dl 

Method 4- Cross- correlation of Instantaneous Phase In- 
crements 

The instantaneous phase (p{t) for both BP and BFV is 
a nonstationary signal and can be decomposed into two 
parts: a linear trend and fluctuations along the trend. 
The trend is mainly driven by the periodic heart rate at 
a frequency ~ 1 Hz. However, the fluctuations are of spe- 
cific interest, since they may be affected by the cerebral 
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autoregulation. To remove the trend, we consider the in- 
crements in the consecutive values of the instantaneous 
phase, defined as 



Aip{U) = ^{U)-^{U-i), (10) 



where ti and are the times corresponding to two 
successive recordings (in our case we have U — U-i — 
0.02sec). The instantaneous phase increment signal Acp 
is stationary in the sense that it has fixed mean and fixed 
standard deviation, and fluctuates in the range (— 7r,7r]. 
In Figs. Etc), [S{d) we show examples of for healthy 
subjects and in Figs. |3g), IS^h) for post-stroke subjects. 

We then apply a cross-correlation analysis to quantify 
the dynamical relationship between the stationary phase 
increments Ai^s of the BP and BFV signals. For each sub- 
ject, during each physiologic condition we calculate the 
correlation coefficient C(t) vs. the time lag r between 
the BP and BFV signals. To quantitatively distinguish 
the control group and the stroke group, we further ana- 
lyze the characteristics of the correlation function C(r). 
Specifically, we investigate the maximum value of C(r), 
denoted as Cmax, which represents the strength of the 
cross-correlation between the instantaneous phases of the 
BP and BFV signals. Another important characteristic 
of the cross-correlation function is how fast the correla- 
tion between two signals decreases for increasing values 
of the time lag r. To quantify this aspect, we choose a 
fixed threshold value, r = 0.3, which is the same for all 
subjects and for which we obtain a good separation be- 
tween the control and the stroke group. Since C(t) is a 
periodic-like function of the time lag r with a decreas- 
ing amplitude for increasing t (Fig.|3, we first record all 
maxima of |C(t)| during each heart beat period (~1 sec), 
then we determine the two maxima with largest positive 
and negative time lags t at which the correlation function 
C(r) is above rCmax- The average of the absolute values 
for these two time lags is marked as a characteristic time 
lag To. 

In summary we note that, to avoid problems with the 
nonstationarity when applying the cross-correlation anal- 
ysis between instantaneous phase increments, we have 
performed the following steps: (1) We have considered 
BFV and BP at the steady state for all four physiologic 
conditions, and we have analyzed only short segments 
(« 2.3 minutes) of data during which the heart rate 
remains approximately stable; (2) We have performed 
cross-correlation analysis after pre-processing BFV and 
BP signals by a high-pass filter, thus removing low fre- 
quency trends; (3) The phase increment signals Ai^ are 
stationary at least according to the definition of weak 
stationarity. For the data segments we consider, Ai^j has 
a constant mean value and fluctuates within a fixed range 
(-7r,-f7r]. 



III. RESULTS 
A. Mean Values 

We compare the mean values of all signals for 
both control and stroke groups and for all four phys- 
iologic conditions (baseline supine rest, upright tilt, 
tilt-hyperventilation and CO2 rebreathing) using the 
MANOVA method. Our results are shown in Table ^ 
We see that the mean values of the BFV and CO2 signals 
are significantly different for the four different conditions 
while the BP mean values are similar. For control and 
stroke group comparison, we find that BFVs from the 
left (stroke side) MCA were significantly different, and 
that the mean value of BP for the stroke group is signif- 
icantly higher than that for control group. Furthermore, 
we observe that the cerebral vascular resistance (CVR) 
is significantly higher for the stroke group. 

B. Transfer Function Analysis 

We apply transfer function analysis on the original BP 
and BFV signals under different physiologic conditions in 
both the low frequency (LF) (0.05-0.2 Hz) and the heart 
beat frequency (HBF) (0.7-1.4 Hz) range. Gain, phase 
and coherence are calculated for each subject and for all 
four physiologic conditions fTable lTHl . We use MANOVA 
to compare our results for the two groups and for the four 
conditions. In both frequency ranges, we do not find sig- 
nificant difference in the gain. In the LF range, phase 
$(/) and coherence 7^(/) are significantly different be- 
tween the physiologic conditions, but are not different 
between the groups. In the HBF range, we find that the 
phase $(/) for the MCAL-BFV is significantly different 
between the conditions (p=0.03). The coherence in the 
HBF range shows no significant difference in the physi- 
ologic conditions comparison, however, it is significantly 
higher for the control group. 

C. Cross-correlation Analysis 

We have performed a cross-correlation analysis di- 
rectly between the BP and BFV signals, after first pre- 
processing these signals by a high-pass filter as shown in 
Fig. n Representative examples for the cross-correlation 
function over a broad range of time lags t for one healthy 
and one post-stroke subject during the four physiologic 
conditions (supine, tilt, hyperventilation, CO2 rebreath- 
ing) are shown in Fig. [5] Group statistics and compara- 
tive tests arc included in Table III. 

The BFV data show more random fiuctuations for 
healthy subjects compared to post-stroke subjects (one 
example is shown in Fig. leading to slightly reduced 
cross-correlation amplitude between BP and BFV for 
healthy subjects. However, the difference in the maxi- 
mum cross-correlation value Cmax is not significant (see 
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Frequency 
Range 


Variable 




Supine 


Tilt 


Hyperventilation 


CO2 rebreathing 


Statistics 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


Low 
Frequency 
0.05-0.2HZ 


Gain 
Gain 


Control 
Stroke 


0.94± 0.30 0.97±0.27 
0.84± 0.42 0.72±0.37 


1.06±0.21 0.97 ± 0.21 
0.95±0.27 0.95 ± 0.25 


0.95 ± 0.41 1.02 ± 0.43 
1.09 ± 0.36 1.08 ± 0.38 


1.03 ± 0.15 1.03 ± 0.15 
0.92 ± 0.20 0.95 ± 0.21 


0.38* 0.11* 
0.48* 0.26* 


Coherence 
Coherence 


Control 
Stroke 


0.65± 0.15 0.59±0.21 
0.53± 0.19 0.44±0.24 


0.75±0.08 0.66 ± 0.19 
0.70±0.20 0.66 ± 0.27 


0.58 ± 0.18 0.62 ± 0.19 
0.62 ± 0.17 0.61 ± 0.19 


0.82 ± 0.09 0.84 ± 0.08 
0.69 ± 0.19 0.70 ± 0.20 


0.0006* 0.0004* 
0.07* 0.07+ 


Phase (rad) 
Phase (rad) 


Control 
Stroke 


0.56± 0.27 0.63±0.26 
0.79± 0.27 0.56±0.38 


0.62±0.32 0.63 ± 0.30 
0.58±0.24 0.60 ± 0.25 


0.92 ± 0.33 0.94 ± 0.39 
0.86 ± 0.32 0.94 ± 0.47 


0.63 ± 0.23 0.63 ± 0.20 
0.42 ± 0.25 0.49 ± 0.23 


0.0002* 0.0001* 
0.74* 0.37+ 


Heartbeat 
Frequency 
0.7-1.4HZ 


Gain 
Gain 


Control 
Stroke 


0.92± 0.18 0.93±0.24 
0.77± 0.29 0.82±0.28 


0.91±0.15 0.91 ± 0.21 
0.88±0.27 0.97 ± 0.43 


0.88 ± 0.23 0.94 ± 0.23 
0.97 ± 0.24 0.98 ± 0.28 


0.89 ± 0.16 0.91 ± 0.18 
0.96 ± 0.28 0.97 ± 0.31 


0.50* 0.73* 
0.96* 0.84+ 


Coherence 
Coherence 


Control 
Stroke 


0.76± 0.15 0.75±0.18 
0.64± 0.24 0.63±0.25 


0.75±0.13 0.71 ± 0.18 
0.57±0.23 0.56 ± 0.21 


0.58 ± 0.16 0.58 ± 0.18 
0.55 ± 0.20 0.56 ± 0.21 


0.68 ± 0.23 0.70 ± 0.22 
0.64 ± 0.25 0.63 ± 0.27 


0.15* 0.27* 
0.035+ 0.05* 


Phase (rad) 
Phase (rad) 


Control 
Stroke 


0.32± 0.13 0.27±0.14 
0.35± 0.31 0.31±0.24 


0.38±0.18 0.38 ± 0.20 
0.36±0.14 0.36 ± 0.21 


0.46 ± 0.27 0.46 ± 0.26 
0.47 ± 0.20 0.49 ± 0.28 


0.46 ± 0.26 0.42 ± 0.24 
0.42 ± 0.21 0.42 ± 0.15 


0.16* 0.03* 
0.95* 0.81 + 



* P value between physiologic conditions comparisons 
f P value between groups comparisons 



TABLE II: Gain, coherence and phase in the low frequency (LF, 0.05-0.2 Hz) range and in the heartbeat frequency (HBF, 
0.7-1.4 Hz) range for the control and stroke group during different physiologic conditions. We compare data from BFV in the 
right MCA (BFV-MCAR) in healthy subjects with data from BFV in the normal side MCA in post-stroke subjects (mean ± 
standard deviation values are presented in the left column for each condition). We compare data from BFV in the left MCA 
(BFV-MCAL) in healthy subjects with data from BFV in the stroke side MCA in post-stroke subjects (mean ± standard 
deviation values are presented in the right column for each condition). The p values from 2x4 MANOVA are calculated for 
comparing differences between groups and conditions. 



Table III), the reason being that the general shape of 
the BP and BFV oscillations at each heartbeat is very 
similar for both healthy and post-stroke subjects, and 
in addition, the amplitude of BP and BFV oscillations 
is much larger relative to the small random fluctuations 
on top of these oscillations. In contrast, the shape of 
the oscillations observed in the phase increments of the 
BP and BFV signals at every heartbeat is very different 
for healthy subjects [see Figs. [SJc), [Sfd)] but very simi- 
lar for post-stroke subjects [see Figs. Elg), EJ^h)]. Thus, 
the instantaneous phase cross-correlation between BP 
and BFV for healthy subjects is significantly different 
from post-stroke subjects, as one can see by comparing 
Fig. 121 and Fig. El This difference is also shown in Table 
III. The comparative statistical tests indicate significant 
difference between the control and stroke group based 
on the instantaneous phase cross-correlation parameter 
C'max when we compare both the stroke side and the nor- 
mal side of post-stroke subjects with healthy subjects. 
There is no significant difference in Cmax obtained from 
the direct cross-correlation between BP and BFV signals 
when we compare the normal side of post-stroke subjects 
with healthy subjects (see also results in Sec. lIIIEll . 

We also note that, healthy subjects exhibit higher vari- 
ability in the intervals between consecutive heartbeats, 
leading to higher variability in the intervals between con- 
secutive peaks in the BP and BFV waves (Fig. As a 
result, both the direct cross-correlation and the instanta- 
neous phase cross-correlation between the BP and BFV 
signals decay faster with increasing the time lag r for 
healthy subjects compared to post-stroke subjects. For 



both methods we show that the parameter tq , which char- 
acterizes the decay in the cross-correlations, separates 
equally well the stroke group from the control group (see 
Table III). 



D. Phase Synchronization Analysis 

In Fig. Owe present the results of the phase synchro- 
nization analysis for four healthy and four post-stroke 
subjects. The statistics for all four physiologic condi- 
tions (supine rest, tilt, hyperventilation, CO2 rebreath- 
ing) in our database and over the entire control and 
stroke groups are summarized in Table III. We find a 
statistically significant difference between the control and 
stroke group based on the entropy index p. 

In the phase synchronization method one estimates the 
difference between the instantaneous phases of two sig- 
nals at the same time. Thus, the phase synchronization 
method does not reflect possible time delays between 
the two studied signals, and correspondingly, between 
their instantaneous phases. Such time delays can often 
be observed in coupled physical and physiological sys- 
tems 1^ ■ To demonstrate the effect of such time 
delay in context of the system we study, we consider the 
difference between the instantaneous phases of BP and 
BFV, when the phase of BFV is taken at time ti and the 
phase of BP is taken at time ^2=^1+ t. As we show in 
Fig. 21 the result of the phase synchronization analysis is 
very different when considering instantaneous phase dif- 
ference with a time delay — the histogram of the phase 
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FIG. 6: Cross-correlation function C(r) of the instantaneous phase increment Aip for the BP and BFV signals during four 
physiologic conditions. We find that the cross-correlation function for all healthy subjects exhibits a very distinct type of 
behavior compared to post-stroke subjects. Two typical examples are shown. (Left) A healthy subject: C(t) has a small 
amplitude at r = and is close to zero at time lags r > 5 seconds during all four conditions. (Right) A post-stroke subject: 
C(r) has a much larger amplitude at r = which lasts for lags r up to 20 seconds, indicating a strong coupling between the 
BP and BFV signals, i.e., loss of cerebral autoregulation. 



difference becomes much broader (the entropy index p 
much smaller) compared to the case when no time delay 
is introduced. 

Statistics for the groups and physiologic conditions are 
presented in Table III, and indicate a very different result 
— no statistically significant difference between the con- 
trol and stroke group for the entropy index p estimated 
for an arbitrary chosen time delay of t = 0.4 seconds. 
This is in contrast to the results obtained for the entropy 
index p from the synchronization analysis when one does 
not consider a time delay. Thus, the results of the tradi- 
tional phase synchronization method strongly depend on 
the time at which the instantaneous phases of two sig- 
nals are compared . This motivates our approach 
to investigate the cross-correlation between the instanta- 
neous phases of two coupled systems at different time lag 
T, as presented in Sec. IIII El 



E. Cross-correlation of Instantaneous Phase 
Increments 

We apply the instantaneous phase increments cross- 
correlation analysis to all four conditions and both study 
groups. We find that the patterns of the cross-correlation 
function C (r) of the instantaneous phase increments Aip 
of the BP and BFV signals are very different for the 
stroke group compared to the control group. In gen- 
eral, the cross-correlation function C(r) for the control 
group is characterized by smaller amplitude and faster 
decay (time lag r less than 10 seconds) [Fig. E]. In 
contrast, for post-stroke subjects, the amplitude of the 



cross-correlation function C(r) is much larger and de- 
cays much slower (over time lags larger than 10 seconds) 
[Fig.E]. Thus for those post-stroke subjects the changes 
of the phase of BFV will change in the approximately 
same way with that of BP signals, indicating a strong 
synchronization. 

The correlations at short time scales (less than 10 
seconds) may be partially attributed to the effect of 
heart rate ('--^l sec) and respiration ('^5 seconds) — i.e., 
they reflect the effect of other body regulations (simi- 
lar to "background noise") on both BP and BFV signals. 
When cerebral autoregulation is effective, because of its 
fast-acting mechanism it may also contribute to the 
significantly weaker cross-correlations at short time scales 
we find in healthy subjects compared to post-stroke sub- 
jects (Fig.|BI). However, the correlation due to the above 
mechanisms will decrease very fast when increasing the 
time lag r between BP and BFV signals. Thus, the cor- 
relations at long time scales (>10 seconds) observed for 
post-stroke subjects (Fig.|H|), cannot be attributed to the 
effect of cerebral autoregulation. Instead, the existence 
of such strong and sustained cross-correlations may imply 
that BFV will passively follow the changes of BP, thus in- 
dicating absence of vascular dilation or constriction and 
impaired cerebral autoregulation for the post-stroke sub- 
jects. 

To quantitatively distinguish the control group from 
the stroke group, we study the characteristics of the cor- 
relation function C(r) for all subjects. For each correla- 
tion function C(r), we first find Cmax, the maximal value 
of C(r), which tells the strength of the correlation. Then 
we choose a threshold value, e.g., r = 0.3, search for the 
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Variable 




Supine 


Tilt 


Hyperventilation 


CO2 rebrcathing 


Statistics 






BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


BFV- BFV- 
MCAR/ MCAL/ 
Normal Stroke 
side side 


direct cross-correl. 


Control 
Stroke 


0.91 ± 0.02 0.89 ± 0.06 
0.92 ± 0.03 0.92 ± 0.03 


0.88 ± 0.05 0.84 ± 0.06 
0.88 ± 0.09 0.89 ± 0.06 


0.88 ± 0.13 0.87 ± 0.12 
0.91 ± 0.06 0.90 ± 0.09 


0.91 ± 0.03 0.91 ± 0.03 
0.92 ± 0.04 0.93 ± 0.04 


0.066* 0.069* 
0.33+ 0.024+ 


direct cross-correl. 

To 


Control 
Stroke 


6.5±4.1 6.6±4.1 
11.4 ±6.6 11.2 ±6.9 


3.3 ±1.4 3.5 ±1.3 
6.5 ±5.6 6.4 ±5.6 


3.2 ±1.2 3.2 ±1.3 
6.2 ±5.8 6.2 ±5.7 


4.8 ±3.2 5.0 ±3.4 
9.8 ±7.2 9.6 ±7.2 


0.008* 0.01* 
0.0001+ 0.0004+ 


phase synch, 
index p for r — 


Control 
Stroke 


0.32 ± 0.07 0.31 ± 0.06 
0.35 ± 0.07 0.34 ± 0.07 


0.25 ± 0.07 0.23 ± 0.07 
0.26 ± 0.09 0.27 ± 0.09 


0.21 ± 0.08 0.20 ± 0.08 
0.25 ± 0.09 0.25 ± 0.10 


0.27 ± 0.06 0.26 ± 0.06 
0.32 ± 0.09 0.32 ± 0.09 


< 0.0001* 0.0002* 
0.047+ 0.008 + 


phase synch, 
index p for r — 0.4s 


Control 
Stroke 


0.18 ± 0.05 0.17 ± 0.05 
0.21 ± 0.05 0.21 ± 0.05 


0.11 ± 0.05 0.11 ± 0.06 
0.11 ± 0.05 0.12 ± 0.06 


0.11 ± 0.06 0.11 ± 0.06 
0.11 ± 0.06 0.11 ± 0.06 


0.12 ± 0.05 0.12 ± 0.04 
0.13 ± 0.05 0.13 ± 0.05 


< 0.0001* < 0.0001* 
0.43+ 0.18 + 


phase cross-correl. 

C,„ax 


Control 
Stroke 


0.55±0.17 0.49±0.22 
0.62± 0.12 0.58±0.14 


0.39±0.16 0.33 ±0.17 
0.43±0.21 0.47 ±0.20 


0.47 ± 0.17 0.43 ± 0.16 
0.51 ± 0.19 0.52 ± 0.20 


0.45 ± 0.13 0.41 ± 0.13 
0.57 ± 0.20 0.58 ± 0.21 


0.006* 0.07* 
0.049+ 0.001 + 


phase cross-correl. 
ro 


Control 
Stroke 


2.0 ±2.0 2.0 ±1.7 
6.6 ±6.7 5.9 ±5.3 


2.2 ±1.3 2.5 ±1.5 
3.8 ±3.4 3.6 ±3.5 


2.0 ±0.9 2.0 ±0.9 
5.3 ±5.6 5.7 ±5.9 


2.6 ±1.5 2.6 ±1.5 
5.9 ±5.4 5.8 ±5.6 


0.6* 0.7* 
0.0003+ 0.0004+ 



* P value between physiologic conditions comparisons 
f P value between groups comparisons 



TABLE III: The maximum correlation strength Cmax and the characteristic lag ro obtained from the direct cross-correlation 
of the BP and BFV signals and from the cross-correlation of the instantaneous phase increments in the BP and BFV signals, 
as well as the entropy index p obtained from the synchronization analysis for the control and stroke group during different 
physiologic conditions. We compare data from BFV in the right MCA (BFV-MCAR) in healthy subjects with data from BFV 
in the normal side MCA in post-stroke subjects (mean ± standard deviation values are presented in the left column for each 
condition). We compare data from BFV in the left MCA (BFV-MCAL) in healthy subjects with data from BFV in the stroke 
side MCA in post-stroke subjects (mean ± standard deviation values are presented in the right column for each condition). 
The p values from 2x4 MANOVA are calculated for comparing groups and conditions difference. 



maximum of |C(t)| during each heart beat period along 
both positive and negative lags, and find two points in 
those maxima with largest time lags at which the cor- 
relation are still above rCmax- The average of absolute 
values of these two points gives the characteristic time 
lags Tq. From tq and C^ax for all subjects and during all 
physiologic conditions, one can confirm that the stroke 
group tends to have larger Cmax and longer time lag tq 
compared to those for the control group. 

We apply MANOVA to demonstrate whether tq and 
Cmax are different for healthy and post-stroke subjects. 
The results are shown in Table ITTll We find that during 
tilt and hyperventilation, the control and stroke groups 
are not significantly different (p values > 0.05). In con- 
trast, during supine rest and CO2 re-breathing, the dif- 
ference between the control group and the stroke group 
becomes significant (p values < 0.05). 

To explain the above difference in tq and Cmax dur- 
ing supine rest, we note that post-stroke subjects ex- 
hibit higher BP and CO2 mean values in baseline (see 
Table P). Therefore, the BP-BFV autoregulatory curve 
for post-stroke subjects is shifted to the right, to higher 
BP values, while at the same time the plateau of this 
curve is narrowed due to the higher level of CO2 CO2 
rebreathing increases the level of the CO2 after the pe- 
riod of hyperventilation (hypocapnia), thus further test- 
ing the reactivity of the C A which shows impaired vasodi- 
latatory responses in post-stroke subjects. In contrast, 
vasoconstrictor responses to tilt-up hyperventilation are 
preserved. 



IV. DISCUSSION AND CONCLUSIONS 

In this study we investigate dynamics of cerebral au- 
toregulation from the relationship between the periph- 
eral BP in the finger and the cerebral BFV for a group of 
healthy and post-stroke subjects during the four different 
physiologic conditions of supine, tilt, hyperventilation, 
CO2 rebreathing. The mechanism of cerebral autoregula- 
tion is traditionally accessed through the response of the 
average cerebral BFV to abrupt perturbation in the the 
average BP (e.g., upright tilt from supine position), and 
is typically characterized by the ability of the cerebral 
blood vessels to restore equilibrium and a steady cerebral 
blood flow after such perturbation — a process which is 
known to operate on time scales above several heartbeats. 
In contrast, we focus on the dynamical characteristics of 
the pressure-flow fluctuation around average values. We 
show that in healthy subjects the CA mechanism is ac- 
tive even during the steady equilibrium state. Moreover, 
we find that on top of the BFV waveforms associated 
each heartbeat there are robust fluctuations which are 
reduced in the post-stroke subjects, indicating an active 
cerebral vascular regulation in healthy subjects on time 
scales within a single heartbeat even under steady BP. 

To test for the dynamical patterns in the BP-BFV 
fluctuations we use four different methods, and we com- 
pare the results of the analyses by evaluating the com- 
bined effects of pressure autoregulation (upright tilt) and 
metabolic autoregulation (hyperventilation and CO2 re- 
breathing) in healthy and post-stroke subjects. We flnd 
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that the gain and phase obtained from the traditional 
transfer function analysis do not provide a signifieant 
difference between healthy and post-stroke subjects. In 
contrast, the coherence is significantly different in the 
heartbeat frequency range (0.7-1.4 Hz) when we compare 
both the normal side and the stroke side of post-stroke 
subjects with healthy subjects. 

Further, we find that the amplitude of the direct cross- 
correlation between the BP and BFV signals does not 
separate the control from the stroke group and between 
different physiologic conditions when we compare the 
normal side of post-stroke subjects with healthy subjects. 
However, comparing the stroke side of post-stroke sub- 
jects with healthy subjects, we find a statistically signif- 
icant difference, suggesting that direct cross-correlation 
method is sensitive to detect abnormalities in CA only 
for the stroke side. In addition, we observe a significantly 
faster decay in the BP-BFV cross-correlation function for 
healthy subjects, reflecting higher beat-to-beat variabil- 
ity in the BP and BFV signals, compared to post-stroke 
subjects. 

Since both BP and BFV are driven by the heartbeat, 

we also test the coupling between these two signals ap- 
plying the phase synchronization method. For healthy 
subjects we observe a weaker synchronization between 
BP and BFV characterized by lower value of the syn- 
chronization index p compared to post-stroke subjects, 
indicating that the CA mechanism modulates the cere- 
bral BFV waves so that they are not completely synchro- 
nized with peripheral BP waves driven by the heartbeat. 
We find a significantly stronger BP-BFV synchronization 
in post-stroke subjects associated with loss of cerebral 
autoregulation. 

To probe how the CA mechanism modulates the BFV 
we investigate the dynamical patterns in the instanta- 
neous phase increments of the BP signals Alpbp and 
compare them with the patterns we find in the instan- 
taneous phase increments of the BFV signals Alpbfv- 
Remarkably, we find that for post-stroke subjects AipBP 
and AcpBFV exhibit practically identical patterns in time, 
leading to very high degree of cross-correlation between 



AifiBP and AipBPV- In contrast, for healthy subjects 
AifBFV exhibits robust random fluctuations very dif- 
ferent from the structured oscillatory patterns we find 
in AifiBP- This leads to a significantly reduced cross- 
correlation between A^psp and AipBFV for healthy sub- 
jects compared to post-stroke subjects. Our results in- 
dicate a statistically significant separation between the 
stroke and the control group when we compare both the 
normal and the stroke side in post-stroke subjects with 
healthy subjects. This suggests that our new approach 
based on cross-correlation of the instantaneous phase in- 
crements of BP and BFV is sensitive to detect impair- 
ment of cerebral vascular autoregulation in both hemi- 
spheres for subjects with minor ischemic stroke. 

Our findings of robust fluctuations in AifBPv, which 
do not synchronize with periodic patterns in Aipsp, 
clearly indicate that the mechanism of cerebral autoreg- 
ulation impacts the dynamics not only on scales above 
several heartbeats but is also active within a single heart- 
beat, i.e., much shorter time scales than previous known. 
Since these fluctuations are present in the data after we 
have truncated the a part of BP and BFV signals cor- 
responding to the initial perturbation related to changes 
in the physiologic condition, our results also suggest that 
the cerebral autoregulation plays an important role even 
in the quasi-steady state. 
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